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Abstract. 

The effect of the heating of neutrinos by scattering with electrons and positrons 
and by e~ e + annihilation on nucleosynthesis is calculated for a spherically symmetric 
baryon inhomogeneous model of the universe. The model has a high baryon density 
core and a low density outer region. The heating effect is calculated by solving the 
Boltzmann Transport Equation for the distribution functions of electron and muon/tau 
VO ■ neutrinos. For a range of baryon-to-photon ratio ln(?7io) = [0, 1.5] and r^ = [10 2 , 10 8 ] 

cm the heating effect increases the mass fraction Xi Hc by a range of AXi Ho = [1,2] x 
10~ 4 . The change of the value of X4 Hc appears similiar to the change caused by 
an upward shift in the value of r]iQ. But the change to deuterium is a decrease in 
abundance ratio Y(d)/Y(p) on the order of 10~ 3 , one order less than the decrease due 
to a shift in 7710- 
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I INTRODUCTION 



When discrepencies arise between observations of light isotope abundances and 
^ the predictions of standard ( homogeneous and isotropic ) Big Bang Nucleosynthesis 

(BBN ) models research has turned to BBN models with inhomogeneous baryon 
distributions ( IBBN models ) to attempt to resolve the discrepencies. Articles 
in the 1970's and 1980's looked for IBBN models with an overall baryon density 
that both equalled the critical density and satirised observational constraints on the 
light elements. [1] These IBBN models could not satisfy observational constraints to 
isotope abundances for baryon densities other than densities demanded by standard 
models. In the late 1990's observations of 4 He and deuterium placed conflicting 
constraints on the baryon density in standard models. ( KKS, 1999 ) [2] found 
a range of agreement for the outer fringes of the 2a ranges of the observations 
using an IBBN model code. Most recently ( KS, 2000 ) [3] used an IBBN model to 
bring baryon density constraints from cosmic microwave background observations 
in agreement with constraints from 4 He and deuterium observations, though not 
with 7 Li observations. 



Nucleosynthesis models have become more descriptive as they include smaller 
effects that faster computers can calculate accurately. One effect is the neutrino 
heating effect. Electrons and positrons pass a fraction of their energy to neutrinos 
through annihilation and scattering. That energy transfer can change the rates of 
neutron-proton interconversion, and then the results of nucleosynthesis. Hannestad 
and Madsen [4] derive a means of solving the Boltzmann Transport Equation for the 
neutrino distribution functions, presenting results for a standard BBN model. This 
author will discuss the neutrino heating effect in an IBBN model. This article will 
show how the heating effect slightly alters neutron distribution and nucleosynthesis, 
and show the effect's dependence on IBBN parameters. 

II THE MODEL 

The IBBN code for this article corresponds to a model with a spherically sym- 
metric baryon distribution. The model is divided into a core and 31 inner shells 
with high baryon density, and 32 outer shells with low baryon density. A run starts 
at electromagnetic plasma temperature T = 100 GK. The distance scale r^, the 
radius of the model when T = 100 GK, can be varied. During a run the number 
density n(i, s) of isotope species % in shell s is determined by the equation [5] 

^~ = MS) £ • 1 ^V M + W m ) 



The first term corresponds to nuclear reactions and beta decays within shell s, the 
second term to the expansion of the universe, and the third term to diffusion of 
isotope i between shells. In this model only neutrons diffuse. As T falls neutrons 
diffuse from the high density shells to the low density shells, until the neutrons are 
homogeneously distributed. The weak reactions n + z/ c ^p + e - ,n-|-e + <-»p-|- 
z7 c , and n <-> p + e~ + v c convert protons in the high density shells to neutrons, 
which then diffuse into the low density shells where the weak reactions convert 
them back to protons. Nucleosynthesis occurs earlier in the high density shells, 
depleting neutrons. Neutrons from the low density shells then back diffuse into 
the high density shells until nucleosynthesis incorporates all neutrons into nuclei, 
mostly 4 He nuclei. 

Figure (la) shows the mass fraction Xi Hc of 4 He nuclei produced by the model, 
as a function of rj and n ( in centimeters ). Contour lines correspond to specific 
values of A4 He . Neutron diffusion occurs later as T{ increases. The shape of the 
contour lines is determined by when neutron diffusion occurs compared with the 
weak reaction rates mentioned above and with nucleosynthesis. As r, increases 
the time of neutron diffusion coincides less and less with the time before the weak 
reactions fall out of equilibrium. More protons then remain in the high density 
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FIGURE 1. X mc (la) and 10 4 AXt He (lb) due to neutrino heating. 



shells. When ri = 25000 cm the innermost shells retain all their protons, and 
nucleosynthesis occurs earlier enough that nearly all neutrons back diffuse into the 
high density shells to undergo nucleosynthesis. X4 He increases as r^ increases to 
25000 cm. But for larger r^ neutrons cannot reach the innermost shells before 
nucleosynthesis, leading to decreasing production of 4 He. 4 He production bottoms 
out at Ti = 790000 cm when diffusion occurs at the same time as nucleosynthesis. 
For higher r, = the high density shells and the low density shells act separately 
from each other, with the high density shells producing large amounts of 4 He. 

To calculate the neutrino heating effect, the model solves the Boltzmann Trans- 
port Equation [4] 
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for the neutrino distribution functions fi{pi). The distribution functions are used 
to calculate the neutrino energy densities, and f Ve (p e ) for electron neutrinos is used 
to calculate the neutron-proton conversion reaction rates. Increased f Uc (p c ) makes 
the conversion reactions produce fewer neutrons for nucleosynthesis, but decreased 
electron energy density p e makes those reactions produce more neutrons at the 
same time. Decreased p c decreases temperature T, making nucleosynthesis occur 
earlier, when more neutrons are present. This last effect, the clock effect, ultimately 
determines the increase in X4 He due to neutrino heating. 
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FIGURE 2. In Figure (2a) the dotted lines are the diffusion rates to and from shell 3, and 
the solid ( increase ) and dashed ( decrease ) lines are the changes of the rates due to neutrino 
heating. Figure (2b) shows n(n) and n( 4 Hc) ( dotted lines ) and the changes in these number 
densities due to heating, n^o is a normalizing factor and a is the expansion rate of the universe. 

Ill RESULTS 



The diffusion coefficient D n is calculated from collisions between neutrons and 



electrons/positrons ( D nc ) and neutrons and protons ( D np ). [6] 
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Neutrino heating affects D nc by lowering the electron number density n c and 
lowering the electromagnetic plasma temperature T at a given time. The neutrinos 
and the electrons/positrons are homogeneously distributed throughout this model. 
So D nc is the same for all the shells. 

Figure (2a) shows the neutron diffusion rates into and out of Shell 3, a high 
density inner shell, for the case of rjio = 3.0 and ri = 25000 cm. Figure (2a) also 
shows the change in the rates due to heating. This change is at first positive, and 
then becomes negative. This change corresponds to a shift in the diffusion rate to a 
time earlier by 0.1 %. All other shells also have this earlier time shift. Figure (2b) 
shows the number densities of neutrons ( n n ) and 4 He nuclei ( n4 He ) at around 
T = 0.8 GK, when nucleosynthesis occurs, for shell 3 and the same values of 7710 



and T{. Figure (2b) also shows the changes in n n and n,4 He due to neutrino heating. 
The clock effect mentioned above is seen in the humps in the graphs for the changes 
An n and An4 He . An n and An4 He are larger between the time when neutrons are 
converted to 4 He nuclei with the heating effect on and when neutrons are converted 
to 4 He nuclei without heating. 

Figure (lb) shows 10 4 AX4 Hc , the change due to neutrino heating. For the ranges 
of i] and rj observed AX4 Hc remains within a range of [1.1,2.0] x 10~ 4 . X4 Hc is 
determined by the baryon number densities in the shells, and those densities also 
determine the magnitude of the clock effect. That seems to account for the contour 
lines of AX4 He tracking the contour lines of Xi He itself. The heating effect on Xi Hc 
seems to be similar to the effect one would get by shifting the value of 77 upwards 
by about 1 to 2 %. The earlier shift in neutron diffusion due to heating would have 
the effect of the lines in Figure (2a) being stretched to higher values of T{. So for a 
given value of r* and rj the change AX4 Hc will have a value corresponding to Xi Hc 
at that point in the graph plus a shift corresponding to a lower value of r,. The 
change is bigger for 7* < 25000 cm, when X4 He decreases with decreasing 7^, and is 
smaller for 25000 cm < rj < 790000 cm when Xi He increases with decreasing tv 

Figure (3a) shows the overall abundance ratio Y(d)/Y(p) between deuterium 
and protons. Figure (3b) shows 10 3 Alog[F(d)/F(p)], the change due to neutrino 
heating. The nuclear reaction rates that create and destroy deuterium depend on 
the electromagnetic plasma temperature T and on shell s's baryon energy density 
Pb(s). Through T neutrino heating shifts the nuclear reaction rates to a slightly 
earlier time. So for r, < 10 5 cm less deuterium remains at the end of a run. This 
decrease in Y(d)/Y(p) due to heating turns out to be one order of magnitude 
less than the decrease that would come from an increase in r)w mentioned above. 
For i] w = 3.0 and r t = 25000 cm one has to increase 77 by 1.53 % to increase 
X*He by as much as neutrino heating does. But that upward shift would decrease 
log[Y(d)/Y(p)] by around 9.19 xl0~ 3 , instead of the decrease of 1.12 xlO -3 done by 
heating. For r* > 10 5 cm the low density outer shells produce considerable amounts 
of deuterium. In that range the heating effect increases deuterium production. The 
neutrino heating effect on 4 He can then be distinguished from a shift in r/10 by 
looking at how Y(d)/Y(p) differs from the results without heating. 

IV CONCLUSIONS AND FUTURE RESEARCH 

The neutrino heating effect increases the mass fraction of 4 He by [1.1, 2.0] x 10 -4 
for the observed range of baryon-to-photon ratio r/10 — [1)32] and distance scale 
Ti = [10 2 , 10 8 ] cm. That's the same order of magnitude as in the standard case. 
The increase in the mass fraction resembles the increase that would come from 
increasing r/ 10 by about 1 % to 2 %. The heating effect also decreases deuterium 
production. But the decrease in deuterium is an order of magnitude less than the 
decrease that would come from the increase of r/ 10 mentioned above. Further, the 
effect increases deuterium production for r, > 10 5 cm. 
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FIGURE 3. log[Y(d)/Y(p)] (3a) and 10 3 Alog[Y(d)/r(p)] (3b) due to neutrino heating. 

A larger future article can show the earlier time shift of neutron diffusion for 
values of i] w and r* other than 3.0 and 25000 cm, and for shells other than shell 3 
as well. These graphs would show how this time shift is characteristic of the whole 
observed range of 7710 and Ti. The article can also look at deuterium production for 
Ti lower than and higher than 10 5 cm, to explain in more detail how the change due 
to neutrino heating changes signs. The article can also comment on production of 
other isotopes, such as 7 Li. One can determine if the change due to heating has a 
distance scale dependence similar to the dependence for deuterium. 

One should note that CMB observations have placed parameters on the baryon- 



to-photon ratio ( Vt B h 2 
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where Vt B h 2 = (3.650 ± 0.008) x lO" 3 ^ ) [7]. 



A future article should particularly look at the IBBN code's behavior in regions 
corresponding to these parameters. 



REFERENCES 

1. Malaney, R.A., and Mathews, G.J., Phys. Repts 229, 147 ( 1993 ) has an exhaustive 
list of such articles. 

2. Kainulainen, K., Kurki-Suonio, H., and Sihvola, E., Phys. Rev. D59, 083505 ( 1999 

)• 

3. Kurki-Suonio, H., and Sihvola, E., astro-ph/0011544. 

4. Hannestad, S., and Madsen, J., Phys. Rev. D52, 1764 ( 1995 ). 

5. Mathews, G.J., Meyer, B.S., Alcock, C.R., and Fuller, G.M., Ap J 358, 36 ( 1990 ). 

6. Kurki-Suonio, H., Aufderheide, M.B., Graziani, F., Mathews, G.J., Banerjee, B., 
Chitre, S.M., and Schramm, D.N., Phys. Let. B289, 211 ( 1992 ). 

7. Buries, S., Nollet, K.M., and Turner, M.S., astro-ph/0008495. 



